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AERONAUTIC SYMBOLS 

1. FUNDAMENTAL AND DERIVED UNITS 


w 

Q 

m 

1 

8 

0 

b 

c 

A 

V 

2 
L 
D 
Do 
Dt 

D, 

G 



Symbol 

Metric 

English 

Unit 

Abbrevia- 
' tion 

Unit 

Abbrevia- 

tion 

Length 

Timft 

1 

t 

F 

meter 

second 

m 

s 

kg 

foot (or mile) 

second (or hour) __ _ _ 

ft (or mi) 
sec (or hr) 
lb 

Force 

weight of 1 kilog;ram 

weight of 1 pound 

Power 

Speed 

P 

V 

horsepower (metric) 

r kilometers per hour 

\meters per second 

kph 

mps 

horsepower 

miles per hour. 

feet per second 

mph 

fps 


2. GENERAL SYMBOLS 


W’eiglit=m^ 

Standard acceleration of gravity =9. 80665 m/s^ 
or 32.1740 ft/sec" 

W 

Mass=y 

Moment of inertia =mP. (Indicate a?ds of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


V Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m“^-s^ at 15° C 
and 760 mm; or 0.002378 Ib-ft”^ sec^ 

Specific weight of “standard^' air, 1.2255 kg/m® or 
0.07651 Ib/cu ft - 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, ^ 

True air speed 
Dynamic pressure, \pV^ 




Lift, absolute coefficient Cx=^ 
Drag, absolute coefficient Cd=^ 
Profile drag, absolute coefficient 
Induced drag, absolute coefficient 
Parasite drag, absolute coefficient 


0 


Cross-wind force, absolute coefficient Dc=^ 


iw 

it 

Q 

U 

R 


a 

€ 

ai 

aa 

y 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 


Reynolds number, 



where Z is a linear dimen- 


sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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SUPERSONIC AIRFOIL PROBLEMS 
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SUMMARY 


LIST OF IMPORTANT SYMBOLS 


A surface integral is developed which yields solutions of the 
linearized partial differential equation for supersonic flow. 
These solutions satisfy boundary conditions arising in wing 
theory. Particular applications of this general method are 
made, using acceleration potentials, to flat surfaces and to 
uniformly loaded lifting surfaces. Rectangular and trapezoidal 
plan forms are considered along with triangular forms adapt- 
able to swept-forward and swept-back wings. The case of the 
triangular plan form in sideslip is also included. Emphasis is 
placed on the systematic application of the method to the lifting 
surfaces considered and on the possibility of further application. 

INTRODUCTION 

The increased emphasis on extending theoretical knowledge 
in supersonic wing analysis has led to a systematic investi- 
gation of the various mathematical methods available for 
treating the basic differential equations. In the present 
report advantage has been taken of the direct analogy which 
exists between the linearized partial differential equation 
for supersonic flow in three dimensions and the two-dimen- 
sional wave equation of mathematical physics. As a result 
of this correspondence, solutions which have been given for 
the wave equation are shown to be applicable to the type of 
boundary condition encountered in wing problems. The 
first section of the report is devoted to the development of 
the solution for the potential of the supersonic flow field. 
The application of this expression to a number of examples 
in supersonic lifting-surface theory illustrates the usefulness 
of such a method of attack. In the first of these examples 
the loadings over the given plan forms are assumed to be 
uniform. The results obtained for such cases appear at 
first to be somewhat academic since undesirable twist and 
camber occur over portions of the resultant surfaces. From 
the uniformly loaded surfaces, however, it is possible to 
develop surfaces having arbitrary load distributions. Im- 
posing the condition that the final lifting surface shall be a 
flat plate leads to the solution of an integral equation in 
every case considered. The results obtained, for some of 
the plan forms considered, have been developed elsewhere 
but not always with the unification of method attained here. 
New configurations are also included among the examples 
given. The methods shown are applicable to a large class of 
unsolved problems of immediate interest. 


a 

A 

b 

c 

C'L 

C{6) 

C, 

E{u, k) 
E,E' 
F(u, k) 


U„lh 

K,K' 


L 

Mo 


Til, no, 713 

Vi 


Pu 

P: {X, Y, Z) 


(I 

S 

So 

U, V, w 

V 

Vo 

A^F,Z 
cn {u,k)\ 
dn {u, k) > 
sn {u, k) ) 
a 

r 

A 


local velocity of sound 


aspect ratio 



span of wing 
chord of wing 


lift coefficient 



load distribution function 

constant value of discontinuity in ip over uni- 
formly loaded lifting surface 
incomplete elliptic integral of second kind with 
argument u and modulus k 
complete elliptic integrals of second kind with 
modulus k and -y l— respectively 
incomplete elliptic integral of first kind with 
argument u and modulus k 
functions introduced in equations (89) and (90) 
complete elliptic integrals of first kind with 
modulus k and ^1— F, respectively 
lift of wing 

free-stream Mach number 
direction cosines of normals to surface S 
static pressure on lower side of lifting surface 
static pressure on upper side of lifting surface 
point at which value of H is to be determined 


free-stream dynamic pressure 



surface enclosing volume V 
area of wing 

perturbation velocities in direction of .V, 1^, 
and Z axes, respectively 
volume 

free-stream velocity 
Cartesian coordinates 

transformed coordinates (See equation (3).) 


Jacobi’s elliptic functions of argument u and 
modulus k 

angle of attack, radians 

Mach forecone from point P: {X, Y, Z) 
semi vert ex angle of triangular wing 


1 
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Ap 

5 

V 

d 

e(u) 

K 

X 

A 

M 


V], V‘2, 
IK?/, 7) 


Po 


( 7 , U 


T 

<Pu 


<Pl 


0 ) 


pressure diirerential (jh—pu) 
angle measured from >V axis 
eonical flow coordinate (See equations (27) 
and (30).) 

/3 tan 5 

Jacol)i’s theta function 

cylinder of infinitc'simal radius enclosing axis 
of forecone V 

surface at whicli strcaun (Mitc'rs induced field of 
wing 

angle of sideslip 

Mach angle of the free stream (^^ = iivc sin ^ 

direction cosines of conormal v to surface S 
incomj)lete elliptic integ]‘al of third kind with 
ai'gument u, j)arameter 7, and modulus k 
density in the free stream 

variables re])resenting either the accelei’alion 
l)otential, the velocity potential, or any of 
the three perturbation velocity components 
surface on which boundary conditions ar(' given 
velocity potential 
acceleration potential 

value of acceleration potential on ui)pei* sid(' 
of lifting surface 

value of acceleration potential on lower side of 
lifting surface 

1-^ 

X 


THEORY 

LINEARIZATION OF DIFFERENTIAL EQUATION FOR COMPRESSIBLE FLOW 


The quasi-linear (i. e., linear in the derivates of highest 
order) differential equation for the velocity-potential ^ in 
the case of compressible fluid flow in three dimensions, is 
expressible in the form 








2^, 




— 




'=0 


( 1 ) 


where a represents the local velocity of sound in the medium 
and Cartesian coordinates are used. Under the assump- 
tions of small perturbation theory (references 1 and 2), 
this e({uation is modified so that it is linear in form and 
consequently more amenable to mathematical analysis. 
Denoting by the variable 0 either the acceleration j)otential, 
the velocity potential, or any of the three perturbation 
velocity components, the linearized expression for equation 

(1) is 

( 2 ) 


where Mq is the Mach number of the free stream and thus 
equal to the ratio of free-stream velocity and tlu' correspond- 
ing speed of sound. 


By means of the affine transformation 

Y= ^±{l-M^)y - 
Z=^±(\-M^z^ 


( 3 ) 


equation (2) can lie put into standard forms. Thus, when 
i\/o< 1 the plus signs are chosen in the radicals of equation (3) 
and e(|uation (2) becomes 

+ = (4) 

while for T/o^l the minus signs are used and, as a conse- 
qiK'iice, 

^xx — ^yy — = 0 (5) 

For the case of subsonic flow (Mq<1) IIr? linearized equation 
is thereby reduced to the well-known Laplace equation in 
three dimensions. Similarly, in supersonic flow (Jfo!>l) 
ecpiation (2) is again reduced to classical type with the re- 
])lacemeiit of the si)ace coordinate A" by a time variable T 
to give the two-dimensional wave equation of mathematical 
j)hysics. The linearization of the general differential equa- 
tion for compressible fluid flow therefore makes available, in 
both subsonic and su])ersonic studies, the results of the ex- 
t(‘iisive work carried out in ])revious research on jiroblems 
related to equations (4) and (5). 

APPLICATION OF GREEN’S THEOREM TO LINEARIZED COMPRESSIBLE 
FLOW EQUATION 

Methods of solution for partial differential equations of 
the type considered here may be classified into two pi'incipal 
categories: methods which express the solutions in terms of 
orthogonal functions and methods which are l)ased on the 
use of (ireeiFs theorem. Volterra’s solution, discussions of 
which may be found in references 3, 4, and 5, applies the 
latter approach to the two-dimensional wave equation and, 
as a consecjuence, his results may be adapted to the study of 
supersonic flow and specific solutions of equation (5). 

Tf the functional notation 

L (12) = 12xx B y y ilzz 

is used, the analytic form of Greenes theorem for equation (5), 
relating a volume integral over the region V to a surface 
ir.t(‘gral over the surface S enclosing U, may be written in 
the form 


J’/J [crZ(12) — 9,L((j)]dV= — J J {(rDn9 — 9Dncr)(lS 


where a, 9 ai’e any two functions which, together with their 
first and second derivatives, are finite and single valued 
throughout the region considered, and 


Dn9 = 7li 


^ ^ 

dX 


where Ui, n 2 , ??3 are direction cosines of inward normals to the 
surface S. 

The expression for Dn9 is, of course, a directional deriva- 
tive. The corresponding term appearing in Green's theorem 
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for Laplace’s differential equation (incoinpressil)le fluid flow) 
is precisely the directional derivative along the normal to the 
surface S. The analogy between the two terms prompts the 
introduction of the so-called conorrnal to S with direction 
cosines ri, ^ 2 , rs defined as 

Vl= Til, V2 = 'tl2, ^3 = 77/3 

The geometrical connection between the normal and the co- 
normal is indicated in figure 1; the angles between the lines 
and the Y and Z axes remain respectively equal, while the 
angles between the lines and the .Y axis are supplementary. 





Kkjukk 1.— goomctric ivlatioiis betwirn normal and conorinal to surfaci* N. 


It follows, in particular, that if the surface S is the .VI" plane 
the two lines are coincident; if S is a cone with semivertex 
angle equal to 45° and axis parallel to the A" axis, the co- 
normal at any point lies long the surface S. 

It is now possible to write 


— Dn^=V\ 


dQ. d9. dQ_d9 


( 6 ) 


VOLTERKA’S METHOD FOR TWO-DIMENSIONAL WAVE EQUATION 


Consider now a surface r which, for the purposes of this 
report, may be thought of as being coincident with the XY 
plane and parallel to the air flow which is in the direction 
of the positive X axis. Two such surfaces are represented 
by the darkened areas in figures 2(a) and 2(b). It is desired 
to determine the value of 9 at the point P:{X^ Y, Z) from a 
knowledge of the boundary conditions given on r. The 
solution to such a problem is immediately suggested by 
equation (8) since that equation requires only the knowledge 

d9 , 

of 9 and — along a surface enclosing a given volume, together 

Or 

with the knowledge of some particular solution a to the wave 
equation valid everywhere within the enclosed volume. 
Further, it is physically evident that contributions to the 
value of 9 at P can come only from points within the forecone 
with vertex at P and also within the envelope of the after- 
cones with vertices at the foremost disturbance points of r. 
Referring to figure 2(a), this would mean the volume bounded 
by the forecone Y and the wedge X springing from the leading 
edge of t; and in figure 2(b), the volume bounded by the 
forecone T and the aftercone X with vertex at the apex of 
the surface r. Since for the boundary-value problems in- 
volved the surface r remains in the plane, equation (8) 
must be applied to all three surfaces X, F, and r. 


Figure 2.— Mach forecone from point P (X, Y, Z) intersecting surface r. 

plan form. 


(a) Rectangular 



XY Plane 


and the surface-volume relation becomes 

[ [ [ [<rL(0)-0L(<r)]rfU= J [[<r dS (7) 

If 9 and a are chosen so as to satisfy equation (5) through- 
out the region V, then equation (7) reduces to the form 

J/" <»> 

The form of equation (8) is a direct analogue to results ob- 
tainable for functions satisfying Laplace’s ecpiation. (See, 
e. g., reference 6, p. 46.) The use of the conormal produces 
this symbolic equivalence. 



Figure 2.— Concluded. (6) Triangular plan form. 
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dil 

Since there is no way of determining Q and along F the 

attempted solution will he especially difficult unless the 
particular solution a and its derivative with respect to the 
conormal vanish everywhere on F. i^ut this is in fact the 
essential part of Volterra’s method of solution. Thus the 
proper choice of <r is 

.Y-X, 


cr=arc cosh 


V(F-T,)^(Z-Z0^ 


(This relation, incidentally, is the indefinite integral of the 
fundamental solution re])resenting a supersonic source in 
three dimensions [(A"— (}"— Fi)^— (Z— Zi)^]~^^^) The 
value of (T is equal to zero on the forecone F since the equa- 
tion of this cone is 

{X-x,y-{Y-Y,y-{z-z,y={) 

and further, since the conormal is always directed along the 
forecone, is the gradient of g along F and is also zero. 

Equation (8) provides an equality for the distribution of 

12 and over X and r, provided 12 and g satisfy equation (5) 

throughout the enclosed volume mentioned. However, al- 
though G satisfies equation (5) everywhere in the enclosed 
volume opposite r from P (under the A^F ])lane in fig. 2), 
along the line (}"— (Z— Zi)^=() (above the AA' plane 
in fig 2) G is infinite and does not satisfy the assunijitions 
made in establishing Green’s theorem. If this line is ex- 
cluded, however, by means of a cylinder k of radius e, with 
axis lying along the line ( )"— Fi)^+ (Z— Zi)^=0, tlien 
equation (8) may he applied to the region outside k and yet 
within the space hounded by X, r, and F. In fact equation 
(8) can then be written 




(9) 


where ti is the portion of r hounding the region of integration. 
If R~ ^ (F— Fi)^+ (Z— Z])^ and cylindrical coordinates e, 
and (A"— A’^i) are used, an element of area on the cylinder k is 
(iS= — ed\l/d{X— Xi)j while 


5(7 5(7 


(A-X,) 


so that 


57? 

-X,)dipd{X-X,) 


lim I' arc cosh clHiX-X,)^ 

f (in •'j--'-’-) ,dH{X-X.) 

--2x (^12(f,F,Z)(7f 

J Xx 


If this result is applied to equation (9), one gets 
and, after differentiating equation (11) with respect to A^, 

aXX X«(“ 

PROCKOURE FOR LIFTING SURFACES AND SYMMETRIC WINGS 

XYhen the region considered is that bounded by the surface 
r, F, and X', the portion of X on the opposite side of t from 
the point P, then g is finite throughout the region and, as a 
direct consequence of equation (9), 




(13) 


where 12' is the value of the potential function on the side of 
T opposite P and p' is in the ojiposite direction to v on r. 
Adding equations (12) and (13), 


“(-V. i',Z)= 2, i, 




1 5 

27t 5A 


i X, s' sXXXX“ 




The integrations over r are now in a form which may be 
interpreted directly in terms of known conditions over bodies 
with given load or symmetrical section. The integration 
over X and X' can be disposed of by discussing the two cases 
shown in figure 2. When 12 is identified with the velocity 
potential, its value can be shown to be zero on X and X' 
regardless of whether the leading edge is swept ahead of or 
behind the Mach cone. When 12 represents acceleration 
])otential or any of the perturbation velocity components, a 
discontinuity exists in the value of 12 for leading edges swept 
ahead of the Mach cone as in figure 2 (a). Analysis of this 
case, however, reveals that for all wing problems the inte- 
gration over X just cancels the integration over X'. When 
the holding edge is swept behind the Mach cone as in figure 
2 (b) the value of 12 is again zero. Thus in any case there 
results the fundamental equation: 


«(A-,y,z)-,Jj,^fXg+“;).rfs+ 


XaXX/“-“' 


) V 

' Op 


(14) 


The counterpart of equation (14) for incompressible fluid flow 
is well known. (See, e. g., p. 60 reference 6.) 

Under the particular conditions for which 


5t2 

dp 


5j2' 

’5/ 


(10) 


(15) 
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O 


over the surface r equation (14) becomes 

The restrictions imposed in equation (15) can be given 
physical significance after the functions 12', and the surface 
T have been given specific meanings. Consider first the 
case where r is a lifting surface. Obviously the normal 
induced velocity w is a, continuous function across r. If 
12 and 12' are velocity potentials associated with the lifting 
surface, 

Y, Z) =w'{X, Y, Z) -g; 

and equation (15) is satisfied. If 12 denotes acceleration po- 
tential or perturbation velocity u, it is necessary to show that 
on the lifting surface 

du_ du' 

This relation holds, however, for since ic(X, Y,Z)=w' (X,Y,Z) 
along r, it follows that 

dw_ dw' 
dZ“dX' 

and from the condition of irrotationality it is possible to ex- 
press the gradient of w in the X direction as the gradient of 
u normal to the surface, that is, in the directions of v and v\ 
Equation (16) is thus applicable directly to lifting-surface 
theory in conjunction with either velocity or acceleration 
potentials. Application can also be made to the determina- 
tion of pressure distribution over the surface of a symmetric 
airfoil at zero angle of attack. In this so-called nonlifting 
case the function 12 is set equal to the induced velocity w, 
T is the plane of symmetry of the airfoil, and equation (16) 
can be used to establish the boundary conditions, provided 
equation (15) is satisfied. For this to be so bwjbv must 
equal —bw'jbv^ But conditions of symmetry give w(Z) = 
— w'{—Z) from which the equality is seen to hold. 

RETRANSFORMATION OF COORDINATES 

Since 

direct substitution into equation (16) yields 

Q{X,Y,Z) = 


n{x,y,z) = 

1 d r r {il—U'){x—Xi)(z—Zi)dxid^ 

2ir dxj J,, [(y— 2/,)H(2— ^i)^— Wl/-2/i)^+(2-Zi)^] 


APPLICATIONS 

GENERAL REMARKS 

Applications in lifting-surface theory may proceed along 
two possible lines depending upon the boundary conditions 
specified. In what is usually referred to as the direct 
problem, or problem of the first kind, the loading is given 
over the wing and the potential function of the flow field 
field is calculated. From the potential function the shape 
of the aerodynamic surface supporting this load can be 
found relatively easily. The inverse problem, or problem 
of the second kind, concerns itself with the determination of 
the loading over a wing surface from a knowledge of the 
surface shape. In the following sections both of these cases 
will be considered. The direct problem will be discussed 
for various plan forms, the analysis proceeding directly from 
the expression for the potential function given in equation 
(17). The detailed discussion of the direct problem is justi- 
fied by its application to the inverse problem where the load- 
ing over flat plates with rectangular, trapezoidal, and tri- 
angular plan forms is determined. The mathematics of the 
inverse problem is less straightforward since the analysis 
involves the introduction of elemental lifting surfaces with 
constant loading and the solution of an integral equation 
for each plan form. 

UNIFORMLY LOADED LIFTING SURFACES IN SUPERSONIC FLOW 

Infinite span wing. — In order to determine the induced 
velocities on the surface of an infinite span, uniformly loaded, 
supersonic lifting surface by means of the methods derived 
in the preceding section, it is convenient to set 12 equal to 
the acceleration potential (p (reference 2) . The lifting surface 
is, in this case, a surface of discontinuity for the function (p 
and corresponds to the surface ti in equation (17). The 
discontinuity in the value of <p between the upper and lower 
surface is equal to 

{<Pu—<Pl)=^ iPl — Pu) 

Po 

where 

Po density in the free stream 
j)i static pressure on lower surface 
pu static pressure on upper surface 


1 f f (12-12') (^-Xi) (Z-Z,)dX,dY, 

27 t bzj J,, [(7- Y,f+{Z-Z,n ^(X~X,f-{Y- Yif-{Z-Z,y 

This solution applies to equation (5) and, in order to relate 
problems to the linearized equation (2), it is necessary to 
use the transformation of equations (3). If the point Xi, 
Yu transforms to the point Xi, Vu if follows that 


It follows that for the uniformly loaded wing in the plane 
01=0 the discontinuity in the acceleration potential is a 
constant, say Co. From equation (17) 


, , Co ^ 



{x—xi)zdxidyi 


( 18 ) 
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Mo 



Figure 3. — Regions of integration for infinite span iinswept wing. 

A sketch of the airfoil plan form is given in figure 3 and 
two possible regions of integration are indicated. In all 
cases the integration with respect to 1 / is performed between 
the limits at which the radical 


^ (i- - (y -]hY+ z'\ 


vanishes while the integration with respect to x depends upon 
the manner in whi(‘h the forecone of the ])oint P intersects 
the discontinuity surface. Denoting the chord length of the 
airfoil by c, the following relations are obtained; 

ip = {) when 

^=±\ Co when i)^xTfiz^c (19) 

ip = () when c<^x^pz 


(When double signs are used, the u])per sign refers always 
to th(' case where and the lower sign corresponds to 

0<O.) 

Th(^ value of the acceleration ])oteiitial is thus seen to he 
zero at all points in s])ace except for those points lying within 
the region between the w('dg('s extending hack from the 
leading and trailing edges of the airfoil. 

It is now possible to determine the induced velocities 
associated with the acceleration potential just obtained. 
Since in linear perturbation theory (rc'fenuice 2), 


dx~ ^ ^ dx dy~ ^ dx’ dz~ ^ ^ dx 


(20) 


where ?/, v, w are respectively the x, ?/, z components of the 
perturbation velocities, it follows that 


1 

y 0 


Tlio induced velocities for the infinite span airfoil result 
immediately from equations (19) and (21). If the upper 
sign of a double sign is again referred to the z'^0 case, the 
r('sults may he written in the form 


Co 

2Vo 


u=± 

v^i) 




>-for 0^x^pz<c 


( 22 ) 


Since the vertical induced velocities are constant, it follows 
that the supersonic airfoil of infinite aspect ratio and uni- 
form load distribution is a flat plate. The relations lietween 
this loading and angle of attack will he considered later. 

Lifting surface with rectangular plan form. —The complete 
discussion of the supersonic lifting surface with uniform 
loading and rectangular plan form is lengthened considerably 
by the fact that in calculating the acceleration potential at 
the ])oint P with cooixiinates x, y, z it is necessary to distin- 
guish between several regions in s])ace in which the point 
may he located. These regions arise from consideration of 
the manner in which the forecone of the point P cuts the 
surface of discontinuity. The value of ip can he found with 
aiiproximately eciual facility in each of these regions hut, 
since this jiaper is concerned ])rimarily with etVects on the 
surface of the airfoil, the solutions for pertinent regions only 
will he given here. 

Figure 4 shows the rectangular plan forrnZA' T' T together 
with the coordinate system to he used. The dimensions of 
the wing are chosen so that the Mach cones extending hack 
from the leading edge will not intersect within the boundaries 
of the wing. This restriction, which is not necessary hut 
merely simjdifies the analysis, imjdies that if h is the span 
of the wing and c the chord length, then 

IV- 1 


\vh(M*(‘ = sin is the so-called Mach angle of the 

Alo 

stream and equal to the semivertex angles of the Mach cones. 
The loading over the loctangular ])lan form is to he uniform 

so the ex])ression ipu—<Pi ^^^'t (‘(pial to f 0 fn*’ ~ "9 ^ 

and oKxiKc. The acceleration ])otential, ('X])r('ssed as a 
function of x, ?/, 2, is thus obtainable from ('(piation (17) and 
the limits of integration must he determiiK'd from the position 
of 1\ From reasons of symmetry, only the ])ortion of sjiace 
for which ?/>0 need he considered. Once tlie acceleration 
|)otential has been calculated, equations (21) may he used 
to calculate induced velocities. The results of such calcula- 
tions arc given and the same convention for double signs is 
used. 

Region 1^: Behind the leading-edge wedge, ahead of the 

( i*ailing-edge wedge, and hounded laterally by the y=0 
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plane and the Mach cone from the leading-edge tip. The 
results in this region correspond to results obtained for the 
infinite span airfoil. Thus 

± 2 (^0 

U= db C0/2T0 

(24) 

Region I 2 : Within the Mach cone from the leading-edge 

tip, outside the Mach cone from the trailing-edge tip, and 
forward of the trailing-edge wedge. Denoting the integrand 
in equation (17) by the symbol /, the expression for when 


y<2 is 

Ax.y. z)~+l (C*. Jr?"''') 


(25) 


where 


^ Xi=x—p (») + 2 ^ 

L=2/+^ ■\l{x—Xi)^—0^z^> X2=xTfiz 

Application of equations (21), after integrating either equa- 
tion (25) or its companion expression when 2/^2 the 

results: 






(y-\ 0 




__Co V 


±2 — &rc tan 






Co 


— Z 






— ~-h/3 arc tan 




[((/-I ») +4-] 


0 ~l 0 




As a partial check of the expression for (p in equations (26), 
it can be seen that in the limit as 2 approaches zero the value 
of <p agrees with the result given in equation (24) on the 
wing while the value is zero off the wing. 

The values of vertical induced velocity in the plane 2=0 
are of particular interest since from a loiowledge of the 
distribution of w the surface shape and local angle of attack 
corresponding to the imposed load distribution can be 
determined. The expressions for w for uniform loading 
will be particularly useful later when the load distribution is 
modified in order to obtain airfoils with specified induced 
velocities. Introducing the notation 




/3 


(27) 


the following results are obtained for the area covered by 
the tip cone: 


For 0<(rj<Cl, 


(28) 


— Co/3 

^^=<'~2xFoJ. V:^ 

and for — (x<^c) 

27tFo V + J_, 


After integration of these two expressions, the explicit value 
of vertical induced velocity throughout the entire region is 
found to be 


w, 


Cop { I Vl—"^^ I • ^ 

-"- 2 ?F.(- 2 + , + "'™V 


(29) 


Equations (28) and (29) indicate that the flow over the 
tip portion of the airfoil is of the type referred to as “conical 
flow.” For this type of flow the values of induced downwash, 
aerodynamic loading, etc., are functions merely of the angle r\. 


831184—49 2 
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Busemanii (referenco 7), Stewart (reference cS), and Lager- 
strom (reference 9) have developed analyses for certain j)lan 
forms which are postulated on the existence of this type of 
solution. In all cases for which the flow field is conical the 
problem is elfectively two-dimensional. Such a simplifica- 
tion reduces the analysis in this report to a consideration of 
a single integral equation while in the references just men- 
tioned complex variable theory can be applied directly. 

Tip of swept-forward lifting surface. -Consider the tip of 
a swept-forward supersonic lifting surface with unifoj-in 
loading (fig. 5), the angle do hetween the leading edge and the 
X axis and the angle 5i between the trailing edge and the x 
axis both being less than the free-stream Mach angle ix. 


Mo 



Figure 5,— Tip of swept forward lifting surface with traces of Mach cones, coordinate system, 
and regions defined for equations (32) and (34). 

In carrying out the integrations it is necessary to distinguish 
between the type in which the tip boundaries are behind 
the Mach cones and the type in which the tip boundaries are 
ahead. The analyses of these two cases are of equivalent 
complexity, however, and can be handled with equal facility" 
by the methods outlined. For all surfaces whose leading 
edges form an apex, only the case where the wing boundaries 
are behind the Mach cone will be considered. A Cartesian 
coordinate system is chosen as shown so that the origin 


lies at the apex, the positive x axis extending downstream, 
the y axis extending laterally, and the z axis being directed 
normal to the plane of the plan form and to the free-stream 
direction. The equations of the sides of the lifting surface 
are 

y = 0 

i p. ^0 

y= —X tan x 

and 

0 

y= — {x—c) tan di=—-^ {x—c) 

The calculation of <p(x, y, z) again must be divided into 
cases de])cnding u])on the location of the point P: {x, y, z). 
In the results listed below are included the explicit expressions 
for <p (x, y, z); the induced velocities, however, are given 
only in the plane z=0, as the integration to obtain a general 
expression is difficult. The velocities in the z=0 plane, 
which are sufficient for the purpose of this investigation, 
can be obtained from a simpler integration since, for tho 
integral involved, 

lim p 

2 — >0 I{x,y^z)dx= 

This simplification was used in the analysis of most of the 
lifting surfaces investigated. As before, it is assumed that 
the discontinuity in <p is equal to Co- Moreover, the ex- 
pressions for Wj^^o Ri’c given in terms of the variables rj 
and 0 ) where 

7,=^Jan<l co=l-^ (30) 


J l{x,y,0)dx 


In this maniKu* the solution is shown to be conical in the 
r(‘gion ahead of the trailing-tip Mach cone (fig. 5). For 
points behind this Mach cone the flow is not conical but a 
function of both rj and co. 

Region /i: Inside the leading-tip Mach cone and ahead of 
the trailing-ti]) Mach cone. Integration of equation (17) 
yields the result 


2tt 


-arc tan 


xy 


w. 


z^,/x^— 


r calculation, 

_ CoISOq 


27tFo, 

1-1 i;i^(’7i + ^o) 

Cofidf) 

f’ 

27tFo^ 



arc tan 


(31) 


(Ivt, 0<7?<1 


Uoi8 rvi— »?' 1 1 1 I 

9 T/ ii cosh I I “b 

SttFoL V Oo. \v\ 


cosh|^^], -!<•,<! (32) 
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Region I 2 '. Inside both tip Maeli cones. The solution in 
this region is simplified through use of the fact that for linear 
diflerential equations any algebraic sum of solutions will be 
another solution of the equation. Since the differential 
e(|uation for the acceleration potential is linear, this property 
can be applied to obtain a solution for the region I 2 by sub- 
tracting from the expressions given for region /i correspond- 
ing expressions in which the variable 5i replaces 8q and (x-c) 
replaces x. Thus 




27T 


— arc tan 






- arc tan 


arc tan- 


(x—c)y 


-arc tan 


xy+x^ ~—^z% 
z^x^—^\y^+z^) 
{x—c)y+{x—cy~Pz-6i 


z^i{x—cY—e^^{y^+^^ 

and for — 1 <t)<1 

_ Cq/ 3 r -y/l— 1_ j 1 I 

“"27rFoL^;~ 00 hr 


z^{x—cY—^‘^{y^-\-z^) _ 
(33) 


2 = 0“ 


- ^ arc cosh ^ - + A arc cosh -i — 

00 i(0o+^)l V 01 1^1 


arc cosh 

01 


* ^4~0i^ I n 
1 0it^+^ 1 J 


(34) 


Lifting surface with trapezoidal plan form. — The linear 
property of the differential equation may be used to advan- 
tage in determining the flow about a trapezoidal lifting sur- 
face with uniform lift distribution, since the boundary 
conditions within the plan form of the airfoil are obviously 
satisfied when the acceleration potential for a triangular 
tip is subtracted from the potential for the rectangular 
surface. 

Suppose (fig. 6 ) the angle of rake of the trapezoid is 8q and 
that 8q is less than the Mach angle /x. The acceleration 
potential will be identical, over the central portion of the 
surface, to that for the lifting surface of infinite aspect ratio. 
Over the parts of the surface which are blanketed by the tip 
Mach cones the flow wull, however, be modified. Because 
of symmetry the determination of this modification need 
only be carried out on one side of the figure. 

If the coordinate axes are chosen as sbown in figure 6 , the 
lateral boundary of the lifting surface is 

y=—x tan 5o= — -- 

It has been shown that both the rectangular plan form and the 
triangular plan form experience conical-type flow over the 
region within the tip Mach cones. Thus, the variable 77 
defined in equation (30) may be used. 

Region /i: Inside the Mach cone originating at the leading- 
edge tip, outside the Mach cone from the trailing-edge tip, 
forward of the trailing-edge wedge, and to the left of the 
y=0 plane. 


For -l<r 7 <(): 


WZ: 


Vi ^1 + 00/ 


J-i 


~iv. [-f +“ ftp 


00 |( 0 o+^)U 


(35) 


Mn 



FiouRE (). — Trapezoidal lifting surface with traces of Mach cones, coordinate system, and 
regions defined for equation (35). 


Swept-back lifting surface. — ^As another example of the 
way in which the linearity of the differential equation may be 
utilized to obtain further solutions, the induced vertical 
velocities for a swept-back wing will be determined for the 
(‘ase in which the leading and trailing edges lie behind their 
respective Mach cones (fig. 7). The boundaries of the plan 
form are given by the equations 

00 00 

2/=-^' (x-c), y=~ (x-c) 

The flow will be conical ahead of the trailing-edge Mach 
cone where the induced velocities can be expressed in terms 
of the variable v- Behind the trailing-edge Mach cone the 
flow will not be conical but will be expressible in terms of the 
c 

variables rj and cj= . 

X 

Consider first the region of conical flow. In order to 
determine Wz=o for a given value of v it is possible to con- 
sider separately the induced effects produced by each half of 
the surface. But in the region ahead of the trailing-edge 
Mach cone, the induced velocities arising from one half of 
the surface are given by the formula for a similar region on 
the swept-forward surface. For reasons of symmetry the 
results for the entire swept-liack lifting surface need only be 
given for values of rj within the limits — l<r 7 < 0 . 
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0 



Figure 7.— Swept-back lifting surface with traces of Mach cones, coordinate system, and 
regions defined for equations (36) and (37). 


In the region where the flow is not conical the solution will 
be built up of a combination of solutions obtained from the 
regions of conical flow. 

Region Ii : Inside the leading-edge Mach cone, outside the 
trailing-edge Mach cone, and to the left of the i/=0 plane. 
For — 1 <^t;<^0 


one of technifiue in integration. The solutions for the uni- 
formly loaded surfaces, however, are particularly useful. By 
methods of superposition these solutions can be used to ob- 
tain the surface loading for specified plan forms (the inverse 
problem) as will be illustrated in the following section. 


L0.40 DISTRIBUTIONS ON FLAT-PLATE LIFTING SURFACES IN SUPERSONIC 

FLOW 


Infinite span wing. — Since the vertical induced velocity is 
constant for the supersonic airfoil of infinite aspect ratio 
(equation (22)) and uniform load, it follows that the airfoil 
is a flat plate. This property distinguishes the infinite aspect 
ratio problem from all other plan forms considered, for the 
load distribution must be modified in the latter cases so that 
twist and camber are removed from the wing to obtain a flat 
plate. 

Denoting the angle of attack of the airfoil by a. 


Moreover, 
and, setting 


it follows that 


'W^2 = 0 ^0 /l/Z' 2 1 

Fo ~2Fo^ ™ 

Vl — Pu = Poi<Pu <Pi)=Pq(^0 

|poFo^ ^ 


(38) 


(39) 


P)liiniiiating Cq between equations (38) and (39), 


Ap_ 4a 

VW-1 


(40) 


fiCodo C” 2 


fiCoOo p 

J-i 


Vl Sq—Vi 


2 dvi 


(30) 


QsL r_ 2 

tFoL ^0 


1 I Vi-gp ^ 


arc cosh arc cosh ^ J 


Region I 2 : Inside both Mach cones and to the left of the 
y= 0 plane. The solution in this region can be produced by 
subtracting from the value of Wz^o given for region Ii the 
value of Wt^Q given for the same region except that in the 
latter case 8 q is replaced by 81 and x by (x~c). Thus, 


Co^ r-2 , 


1| , Vl-0o' ^ (2-eo^-7j2) , 

^1+ aiccosh + 


~ arc cosh 


di 


arc cosh 




(2-0^2) 

-d,W) U 


(37) 


Although the uniformly loaded lifting surface was the only 
prescribed loading analyzed, it should be noted that the 
basic integration leading to a solution of this type of ])rob- 
lem (equation (17)) is in no way restricted to a uniform load. 
Arbitrary loadings that may or may not be analytic functions 
of X and y can be specified and the problem therefore becomes 


The result given in equation (40) is the well-known Ackeret 
expression developed in reference 10. The derivation here 
follows the approach of Prandtl (reference 2). 

Rectangular plan form. — Since the vertical induced velocity 
for the uniformly loaded supersonic airfoil of rectangular 
plan form is not constant over the portion of the wing 
covered b}^ the tip Mach cones, it is necessary to modify the 
load distribution within this region in order to get a flat plate. 
The determination of the required load distribution will be 
shown to depend on the solution of an integral equation and 
subsequent problems dealing with other plan forms will, 
from a mathematical standpoint, be similar in form. 

The rectangular plan form will be thought of as being built 
of superimposed trapezoidal lifting surfaces with variable 
angles of rake (fig. 8), each trapezoidal surface having a uni- 
form load distribution but with loading allowed to vary with 
the variable rake angle 8 . 

Since the flow over the part of the airfoil within the Mach 
cone is conical, it is possible to express as a function of 
7 ] where 



Setting 

P tan 8 = 6 
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and using equation (35) 

p-'p p 

P Je=o \ J-i Vi Vi-r^/ 

where C'(6) = (pu—(pi for the single trapezoidal surface 
rake angle 5. 


(41) 

with 


Mo 



Vi 


Figure 9.— Region of integration showing line of singularity for equation (44). 



It can be shown from ecj[uation (35) that, if C'(6) is a con- 
tinuous function, 


lim C"(— 7? — €)(?i(r), — 7) 

L 




Hence 

c'{B)de 

Vl Jo 


(42) 


and the solution of this equation is 


X 

Figure 8.— Rectangular plan form built of superimposed trapezoidal lifting surfaces with 

variable rake. 

The solution of the problem depends on the determination 
of a function C' {6) which, when substituted in equation (41), 
will yield a constant value of Wz=o{tj) ; that is, a value of 
independent of the variable rj. Imposing the condition that 

dw,^o_r. 

drj 

the problem is resolved into one of solving the equation 


0=4 r~'c'{d)dd r ’'-i' 

avJt=o J-i Vl 


171 + ^ 


C'{e)= 




(43) 


where Ci is a constant to be determined later. Substituting 
from equation (43) into equation (41) 


2TrVoW{iri) ^^0 

pa 


irdS 

p de r 

jo v0(l-6>) 



V 1 — 


drji 


(44) 


The region of integration in the rjid plane for the double 
integral of equation (44) is shown as the cross-hatched area 
of figure 9, a singularity in the integrand occurring along the 
line d= — 7]i. Rewriting the equation and reversing the 
order of integration in the double integral. 


By means of the notation 



dvi 


= G,iv,d) 


the integral equation is written in the form 

o=lim |~4 c'(e)Gi{v,d)dd+^ f C' {d)G,{v,e)dd] 

€-40 L^'nJe=o J 


where the singularity in the integrand necessitates the use 
of the infinitesimal e. The evaluation of the derivative 
thus leads to the expression 


C'(d) 




C'{—v — i)Gi(v, — V — G + C'{—V + i)Gi(r], — ri-\-t)~^ 


)=lim r f” 
LJo 


i)Gi 






■ 2x arc sin 


inV^l — lim f — — di)i 

Jo Vl 


r f-'- a + p 

L^o V^(i — ^) J~’ 


dd 


’?+« (^1+^) V^(i 


=] 


The bracketed expression in this equation can be shown to 
vanish for all values of rji between zero and —1 so that, 
finally, 

— CiPtT 


Wz=o ^ 2 t ai 
ZtT Ko 


;’c sin J ^ | = 

Jo 


2Vo 


(45) 


Since the trapezoidal lifting surfaces are superimposed, 
the loading C{6) over the resultant rectangular plan form 
satisfies the relation 


dC{d) 

dd 




a 

V^d-e) 
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Imposing tiio condition that (\d) = 0 at 6—i), it follows tliat 
C{6) = 2C\ arc sin ^fd (46) 


This equation gives the incremental change of acceleration 
potential between the upper and lower lifting surface of the 
rectangular wing. As a result the increment in pressure is 

Vi—Pu = Po{^u—^i) = 2poCi arc sin 


Expressing tlie jiressure difference in nondimensional terms, 




.Pol^V 




arc sin 




(47) 


The constant (\ may be eliminated between equations (45) 
and (47) and as a consequence 


(Z 


W; 


g 

- 7 , arc sin \6 
ho i^TT 


Since the angle of attack a of the airfoil is by definition 
equal to — ^‘V^;the final expression for the loading, in coefli- 
cient form, over the outer portions of the rectangular wing is 


airfoil must be adjusted properly to give constant induced 
vertical velocity. Superposition of trapezoidal lifting sur- 
faces with loadings varying with rake angle 5 can again be 
used and the conical nature of the flow employed. Setting 

rj = ^ylx 

6=^ tan 8 

6q=P tan 8q 


equation (35) leads to the expression 


2jryo^Mi=.o= c>ie) r de 

P Jd=d0 \ J-1 Vi 

(52) 

where C' (d)=(pu—(pi for the single trapezoidal surface with 
rake angle 8. 

The analysis in this case follows along lines directly analo- 
gous to that used for the rectangular surface. For the 
present configuration the loading function for the superim- 
posed trapezoids is given by the relation 


C'(d) = 


a 


(53) 


<z 


8a 

TrV.W-1 


arc sin 


(48) 


and the integration to obtain Wz=o can be simplified to 
give, as a final result. 


The general apiiroach used to obtain this result is similar to 
that used by Schlichting (reference 11). The error in 
Schlichting’s final result has been noted by Busemann and 
others. 

Lift coefficient Cl for an arbitraiy wing is defined by the 
relation 


Cl — 


L 

Soq 



(49) 


/ 3 < 7 i r* irde 

2 ;rFo J.0 y{e-e,) {\-e)~ 2Fo‘ 

The loading C{d) over the resultant trapezoidal plan form 
can be found from the relation 

dC{d) C, 

^ de V(0-0o) {1-6} 


where 

Z = total lift of the wing 

element of area on the wing 
Sq= total area of wing 

For the rectangular wing the values of Ap/q over the tip 
and center sections are given by equations (48) and (40). 
As a result of this integration 



where A is the asjiect ratio and by definition equal to the 
ratio of the siiuare of the span and the wing area. As a 
final conclusion the lift-curve slojie of the wing is 


From the boundary condition that C{d) = 0 at d=do it 
follows that 


and 


C{d) = 2Ci arc sin 


le-Bo 

\ 1-00 


Pl-Pu 
I P,JV 


Ay>_4F, 
2 “Fo^ 


arc sin 


V' 


e-d o 

1-00 


(55) 


Elimination of C\ between e(iuations (54) and (55) and 
introduction of angle of attack a for — Wj^o/Fo gives as 
aerodynamic loading over the portion of the airfoil within 
the tip Mach cones the expression 


dC^^4 

da P 



(51) 


Ap 8a 

g tP 


arc sin 


/0-0O 

Vl-0o’ 


O<0o<l 


(56) 


Trapezoidal plan form. — The results given in ecpiations 
(48) and (51) are capable of generalization to the case of the 
flat plate having trajiezoidal plan form and with rake angle 
8o less than the Mach angle of the stream. For such a 
configuration the airfoil is again blanketed in part by the 
tip Mach cones and the loading in this outer section of the 


Figure 10 indicates the variation of the loading over the tip 
section of the trapezoid. The variable (^/a){Ap/q) is 
plotted against ^ tan 5 for jS tan 6 q equal to 0, 0.3, and 0.6. 
The curve for /3 tan 6 q= 0 corresponds to the case of the 
rectangular wing and shows results in agreement with equa- 
tion (48). 
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Figure 10— Load distribution over tip for various trapedzoidal plan forms. 


By moans of oquation (49) together with equations (56) 
and (40) the lift coefficient of the trapezoidal wing is expres- 
sible in the form 


C^= 


4a 


c c 

^~2b 


1 — ^ tan do 



(57) 


Introducing the aspect ratio A of tlie wing where 

■^ = —7 ^ 

c f 1 —■ j tan So) 

one gets for lift coefficient tlie relation 
From equation (57), 


dC^^4 

da 

.dC,,. 


c c 

1 — tan 60 - 9 ^ tan ^ 


2b 


1 — j tan 5o 


(58) 


(59) 


In figure 1 1, /3 is plotted as a function of Afi for ^o=0, 9’ 



Kir.rRE 11.- Variation of roduced lift-oiirvo slopo with rtnhiood a.^!pcct ratio 8 A for 

da 

various tra[)ozoidal plan forms. 


terminated at values of Aj3 for which the tip Mach cones 
intersect on the trailing edge of the wing. 

Triangular plan form, type 1. — The pressure distribution 
over triangular lifting surfaces with constant induced vertical 
velocities will be developed in the following three sections. 
These plan forms are indicated in figures 12(a), 12(b), and 
12(c) and shall be denoted, respectively, as types 1, 2, and 3. 
Types 1 and 2 are actually special cases of type 3; namely, 
the cases where one leading edge is parallel to the free stream, 
and where both leading edges make equal angles with the 
stream direction. Type 3 includes any plan form which has 
leading edges swept behind the Mach cone but on opposite 
sides of an axis drawn through the vertex of the triangle 
and parallel to the free stream; and, further, has a trailing 
edge such that the Mach cones from either tip do not cross 
the surface of the wing. The principal reason for considering 
the three types separately is to show the manner in which 
the span wise loading appears in the solution of the problem. 
In types 1 and 2 the proper load distribution is found readily 
while the final type requires a more careful treatment. 

In order to determine the load distribution over the airfoil 
it will be convenient to use a differential element over which 
the loading is uniform. The elements may then be summed 
and the distribution of loading adjusted so that the induced 
vertical velocity at any point on the total lifting surface is 


and 1. The curve for agrees with results given by 

equation (51) for the rectangular wing. All curves are 
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Mo 



(b) 


Mo 



(o) 

Figure 12.— Trianpiilar flat platu lifting surfacos. (a) Type 1. (6) Type 2. (c) Type 3. 

constant. For the triangular plan forms it is possible to 
assume that conical flow exists and the analysis may be 
carried out using the angular coordinates that have already 
been introduced. 

Figure 13 shows the elemental lifting surface to be used. 
The sides of the element extend back from the tip of the 
Mach cone; making angles 8 and 5 + A5 with the positive 
X axis or free-stream direction. Corresponding to previous 
notation, the relations 0 = ^ tan 8 and d-\-A6=P tan (5 + A^) 
are used. The vertical velocity induced by the element of 
surface may be denoted by Aw and it follows that 

Aw=w{6-{- Ad, rj)—w{6, rj) 


where imdw{d+A6,r]) are the velocities induced by the 

triangular-tip surfaces with uniform loading and with tip 
angles equal to 8 and 8 + A8, respectively. Applying a 
limiting process, 


lim lim Me + M,r,)-wie,v )^^ 


(60) 


It follows that Wz=Q for the resultant lifting surface will be 
evaluated by an integration with respect to d. If dWz^ojdd 
can be expressed in the form of an integral with respect to 
7 ], the relation for Wz^q will then be similar to those given 
in cfiuations (41) and (52) for the previous plan forms and 
the expectation will be that the function C{d) can be deter- 
mined to give constant induced vertical velocity. 

The method of attack just outlined is postulated on the 
existence of an integral expression for dWz=o/^d. Such an 
(expression is, however, obtainable directly from the integrals 
in eciuation (32). Integrating tluese relations l>y parts 
after first differentiating by d, leads one to the formulas: 


For — 1 <^r]<^0 


§Co r -y /l— , p ^ ' 

dd 27rVo L’?(^“l“^o) J-i ’7i^(’7i + ^o) Vl 

and for 0 <r 7 <l 

/3Co r y_i I p 

bd 2x1^0 Ji ’)i^(’7i + ^o) Vl 


(61) 


(62) 


If the elements are summed over the type 1 triangular 
wing, induced vertical velocity is 




(63) 


where C(d) = <pu—<Pt (oi‘ the element at 5=arc tan If 
Wz=Q is constant, then 

dv 

and from this criterion the function C(d) will be determined. 
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Thus, using methods similar to those introdiieed in the 
develo|)ment of equation (42), 


0^ 


d 

dr] 


V jo V 


Cid) 


dd-\- 2 


77^1- 


vj: 


C{6)de 

n+0 


ro“C{ 
Jo V 


0 “ C{e)de 


(64) 


Tl»c general solution of this e(ination is 


r(0) = 


c,e+c, 

^d{6a-d) 


where Ci and Ci are constants. Since, however, the Kutta- 
Jonkowski condition requires that loading vanish along the 
edge 6=0, it follows that (72 = 0 and the required loading 
takes the form 


C{6) 




e 


(65) 


If equation (65) is substituted into equation (63), vertical 
induced velocity can he calculated from the expression 


2tVqWz 


/3(7, 


V jo 

no ode n 

Jo J-i 


(r] + 6) V^(^o — ^) 


dvi 


+ VI —vd 


(66) 


The region of integration in the rji, 6 plane for the double 
integral of ec[uation (66) is shown in figure 14 for the case in 
which —do<^r}C^O. A singularity in the integrand of the 
doul)le integral exists along the line 6-\-r]i = 0. Reversing 
the order of integration, equation (66) may he rewritten as 


.r: 


de 


de 

{6+r,^)^^ 


(0+7,) 

dvi 

-floT/i Vl 


r~ea 

J-i 


==— f’ - p' 


di)i 

de 


(0+^?i)Ve(0^) 

(07) 


The single integral in e((nation (67) has a singidarity at 
6 = —7, since — 0o<C’7<CO and t, therefore lies inside the region 


e 



Fir.uRK 14. Region of inU'gration showing line of singularity for equations (60) and (78). 


of integration. A corres])onding singularity occurs in the 
second of the double integrals at ^= — 771. Consider, there- 
fore, the integral 


r 


dd 


(^+v)^e{^o- 


— =lini 

- 6) t^fl 


u: 


de 


£ 


(0+ rt) -\j 60 0 

de 


+ 


l7+« (0+7,) V^(^t 




( 68 ) 


The indefinite integral is 


1 — 77 ^q-|- ^^Q- h 2?7^ — 2-\A( — 77^0 77^) ^(^0 ^ 

0+77 


SO that the definite integral is 


lim 

€->() 


1 



do 

-do 


+ 


[— 277^0— €^— 277^—2677 — 2 V(— 77^0— ? 7 ^)(— 77— e)(^o+ ^ + e)] € ) 

[— 277^0 H-e0o—277^-h2e77—2-\ (— 77 ^ 0 — '^^(“'^+0(^0+^— e] (— e) i 


The value of this expression is 0 and eciuation (07) therefore 
becomes 


2ttVo'W,^o_ _ ^^d ^ 

&C\ J-i 771^1 — rjdj^ (d-{- r]i) \ dido 6) 


(69) 


Since, in this region of integi'ation, — 1 <^t 7<^ — it follows 
that 


n dd 

Jo {d+r],)^d{d,-d)~ 

1 — 77i^o4-^^o+ 277i^ 

'*‘”2V(0iU’n0o)e(0o-0)Jo 

and 

t^Ti 

^ 2I 0J-1 771-vl — t7i‘^-\ 77i(t7i-1-0o) 

The integral of equation (70) can be transformed by means 
of purely algebraic substitutions into a form that integrates 
immediately into complete elliptic integrals of the first and 
second kind. However, in the consideration of the type 3 
plan form it will be necessary to resort to other methods of 
transformation, so that a more uniform approach, employing 
Jacobian elliptic functions, will be used throughout. (vSee 
reference 12.) 

The quartic under the radical in equation (70) is first re- 
duced to an expression of the type appearing in elliptic in- 
tegrals of canonical form. This is accomplished by succes- 
sive application of the transformations 

k -|- Is I k 

’6 =-T+s aiids=y 

where k and I are cliosen so as to destroy the odd powers of 
the variable. By means of these transformations, induced 
velocity becomes 


— 7T 

\ 77(^0+^) 

(70) 
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, (l-kW 

2V"o Vit(«o— A-)(l-A:2)Ji 


where E' is tlie complete elliptic integral of the second kind 
with modulus k' = Equation (72) can be further 

simplified by writing k in terms of 6q so that 


where 

i-VT^ 

00 


The integration of equation (71) will be performed after 
first considering two parts such that Wz=q=Wi-\-W 2 where 

'2Fo V^*(0o-^)(1 

and 

_ /3(7, (1— P)F {\—f)dt 

2\\ ^kie.-kyn-k^)}, 


This separation is prompted by the fact that the integral for 
Wi is expressible in terms of elementary functions after the 
simple transformation k^=z. The results of such an integra- 
tion lead to a value that is zero at the lower limit and infinite 
at the upper limit. However, an inspection of the original 
integral in eejuation (70) shows that is finite so the 

infinity obtained for Wi must be canceled liy a corresponding 
infinity of equal magnitude in Wo. The actual proof of this 
statement necessitates, of course, treating the combined 
expressions as an indeterminate form where the upper limits 


of the integrals for Wi and W 2 are replaced by y+ e and the 

rC 


limit is taken as e approaches zero. 

Introduce now in the integration of W 2 Jacobian elliptic 
functions and set 


so that 


k) = snu 


df — enu dnu du 


The expression for W 2 becomes 

_ fiC, (l-P)K cn^u 

2Fo ^lk{do—kJil—Fj Ja dn-u^*^ 

where K and K' are the complete elliptic integrals of the 
first kind with respective moduli k and k'=\\—k^. Inte- 
grating and combining with iVi, one has 


W,^Q = Wi + 


{l-E)kH 

2Fo ^■k{eo-k){i-k2) 



snu enu 
dnu 



R+iK' 

K 


where E{u) is the incomplete elliptic integral of the second 
kind. After substitution of the limits, induced vertical 
velocity is 


^z = 0= — 


pc, {\-k^)E' 

2Fo \k{do-k)(l-k'^) 


w, 


pc, /2(l + Vl-0o^) p, 
2Fo V ^ 0 ^ 


(73) 


where the modulus of E' is ^Jl—k‘^ and k 
For the loading in question 

Vl — Pu = Po((pu — ^l)=PoCi 0^ 




and 


Pi—Pu ^^p^ 2 Ci I e 

\poV^ ^ Vd‘\0o-t 


By means of eciuation (73) the constant C, may be elimmated 
and 


Ap_ 2a /2<9(1 — Vl — V) 
q ~I3E' \ 00-0 


(74) 





(72) 


Figure 1.5.— Load distribution over triangular plan forms of type 1. 


VOLTERRA’S SOLUTION OF WAVE EQUATION AS APPLIED TO THREE-DIMENSIONAL SUPERSONIC AIRFOIL PROBLEMS 17 


Figure 15 shows the variation of ~ ~ with /3 tan 8 for values 

of tan do equal to 0.3, 0.6, and 0.9. 

From equations (74) and (49) the lift coefficient of the right 
triangle wing with trailing edge normal to the free-stream 
direction can be determined. It follows that 


— Vf - 


(75) 


Since the aspect ratio A is given by the equation A=2Bolfi, 
equation (75) can be used to find the lift-curve slope as a 

dC 

function of A. In figure 16 a plot of ^ as a function of 
pA is given. 


or, since C{d) = C{—6), 


y^iv)z=o T/ 

7T V f 


(76) 


The function C{6) in equation (76) must give a constant 
value for ^^;(77)^=o so that will vanish. Imposing 

this condition it can be shown that a solution is given by 
the relation 

Cl 


and, after substitution in equation (76), 


(77) 


Mo 


t 



tVoW^=o_ P® ^il—rj^dd _ 

/3Ci “ Jo 

de p rfr/i 


(78) 


The integration of equation (78) is to be performed under 
the assumption that — 0o<C'^<O so that the region of 
integration in the r]i,d plane for the double integral is as 
shown in figure 14. Reversing the order of integration in 
the double integral, eciuation (78) may be written in the 
form 



Figure 16.— Variation of reduced lift-curve slope fi with reduced aspect ratio /3.1 for 

da 

plan forms of type 1. 


ttVoW,^ 

pc] - 


dvi r 

J-l 7/1 Vl —'^1^ do 

p r 

J-eo r}i^l—7]d Jo 


de 


de 


^0 de 


0 {e^-r^d)^eo^-e^ 
Evaluation of integrals of the form 

po 

is accomplished by means of the substitution 


Triangular plan form, type 2.— Figure 12(b) shows the 
symmetrical type of triangular plan form considered in this 
section. The semivertex angle is 8o and Bo is defined by the 
relation 

Bo=fi tan 5o 

The loading element used in the previous section can be used 
again and equations (61) and (62) are applicable directly. 
Because of the symmetry of the figure, it is necessary merely 
to insure the constancy of Wz=o over the left half of the wing 
in order that the entire wing be a fiat plate. 

Summing the elements over the type 2 triangular wing, 
induced vertical velocity over the portion of the wing for 
which — 1<^77<^0 is 



After substitution, the integral becomes 



and, by straightforward integration, 

0 forf,<l 


1= 


9*, /v,2 f) 2 n ; 


for ^- 2>1 


(79) 
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This result shows that the second double integral of equation 
(79) vanishes (since as does also the single integral 

in the equation. For the remaining double integral, however, 

and 


Setting 


w 




dy)\ 


Vi 


S/1- 


~Vi —^0 


Z = 


1 


(80) 


equation (71) transforms to 


^Ci 

2T^oJi V(^^-l)(l-^oV) 


Introducing the modulus k = 6o and making the substitution 
z=sn{u, k) =snu 


one gets the expression 


^ 1 ^ 1 = 0 = c^rr sn^udu 

oJ K 


ilsa 

2V,k^ 


\u-E{u)] 


K+\K' 

K 


2F„F 


E' 


wliorp the prime again refers to tlie (•om])lementary modulus 
of the complete ellij)tic integral. Since k = 6o 

W^=0=~9pS-2^' (81) 

0^0 

where A:' = Vl — 

For the loading in (piestion 

Vi—Pu_2,p^ 2(\ 

^ nV^ IV A ^ 0 ^ — 

9 Poy 0 


so that, eliminating C\ between this equation and equation 
(81), 

Ap_ 4 qiV 


2 

Ap 


(82) 


Figure 17 shows the variation of — ^ with /3 tan 5 for values 

® a q 

of jS tan do equal to 0.3, 0.6, and 0.9. 

From equations (82) and (49) it is possible to find the 

expression for lift coefricient of a triangular or delta wing. 

Thus: 

2tv6ool 

■ m' 


Cl = ^ 


(83) 


Since aspect ratio of the wing is 


Mo 




/ 1^)3- 

where the modulus of PJ' is k' = -^1 1 - ^ • This result 

agrees with that obtained in another manner by wStewart 

dC 

(reference 8). In figure 18 a plot is given of ^ dJ^ ^ 
function of pA. 

Triangular plan form, type 3. -Figure 12 (c) shows the 
plan form now to be considered. Relative to the x axis or 
free-stream direction the sides of the triangle form the angles 
5o and 5i so that the total vertex angle is 6o + 5i=^2A. The 
variables 6q and 6i are also introduced satisfying the relations 
6q=^ tan 5o, 0i = j3 tan 5i. The same loading element that 
was used for type 1 and type 2 triangles may be used and 
equations (01) and (62) apply. It will then be necessary to 
determine the distribution of load so that the induced vertical 
velocity over the plan form is a constant. Since this induced 
velocity must be the same on both sides of the 5 = 0 axis, 
two equations result: 


.400 

^ tvolA 

^ L-2PJ' 


(84) 


For — 1 <? 7<0 

- 0 = 2 V J_,. [fel + J_,' -2 


(lr)t 


+ a' 1 — 7/lU 


1 C{d)de 


(85) 


lift coefficient becomes 
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3/o 

V 




Figure 18— Variation of reduced lift-curve slope d with reduced aspect ratio for 

(la 

plan forms of type 2. 


and for 0 <^? 7<^1 




27t VoJ-Oi \__r} {rj 6) 



dr] I 

vd{r]i~\~^) A 1 — Vi 



(86) 


Tlie evaluation of ^o) ^^nd IhiPxy ^o) is accomplisliod 

ill the same manner as has been used previously: first, a 
reduction of the quartic under the radical to canonical form 
and second, transformation by means of Jacobian elliptic 
functions Jollowed by direct integration. Since the cacu- 
lations for both equations are quite similar, only in the case 
of ^o) will the details be mentioned. 


a-\-hs 


By means of t lie transformations k]y= . , and where 

' l+s t 


_ l-Mi-Ad-0(l-^i^) 

^0+^1 


A_ l-Mi+A(l-^o^)(l-^d 

^0+^1 

The introduction of the symbols R and k defined as 

h—a 


H= 


V ( 1 — (^ 0 — (k) { 6 1 + a) 


k = 


e„-b 


oquation (89) reduces to 


Ih = amk 


r\ 

Ji (]— 


-hi(a — b) — l]dt 




(91) 

(92) 

(93) 

(94) 

(95) 


Tiio inlegi'and divides iiatui'ally into two parts, one containing 
even powers and one containing odd ])owers of t in the nu- 
merator. The latter j)art integrates into elementary functions 
after substituting f=u^ and equation (95) thereby becomes 


From the solutions to the problems of type 1 and type 2 it is 
possible to construct a solution of the more general problem 
by expressing the loading function in the form 


C(e) = 

V(d,-h0)do-0) 


(87) 


where A and are constants that can be determined in terms 
of Wz=o from equations (85) and (86). Equation (87), in 
conjunction with equations (85) and (86), yields the expres- 
sions 

^i) “b7^/^2(^0) ^i)] 


where 




2Vo 


[-AliA0i,eo)+BH2(e:,do)] 



dvi 

viV(i— vd) (vi— ^i) (^1 + ^o) 

dvi 

7}d ( 1 — vd) ill — ^i) (^1 + ^o) 


( 88 ) 

(89) 

(90) 


Setting 




(f—\)dt 


and substituting 


one gets 


x=sn{u,k)=snu 

I 2 — ^1 l-[“(^ 1) "i 2 T~ 

JK L l—a^sn^u_\ 


du 


If sn h BOW may be written as 

j __ . dny r^+iK' Ic^snycnydnysn^u , 

^ ^ k^snycny Jk 1 — k^srdysrdu ^ 

or 

(97) 

where IT ( 1 ^, 7 ) is the fundamental elliptic integral of the 
third kind. 
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The evaluation of U(u, k) is best acliieved by moans dfjits 
expression in terms of thcla functions ami zota functions. 
Thus 


II K 7) = 2 log 


Q{u—y) 

Q(u+y) 


{-uZ{y) 


and the brac'keted term in equation (88) is 


U{K+iK',y)-n{K,y) = ^ log 


0 {K+ il^ — t) ^ t) I 

0 (K— 7) 0 (K+ iK ' + t) ^ 


iK'Z{y) 


The theta functions are qiiasi-perioclic, that is, they satisfy 
the relations 

Q{u-]-2K) = Q{u) 


From this property, together with the fact that 0(?/) is an 
even function, it follows that 


0(7v + i/^'-7)0(/l + 7) 
0{K-y)e {K+ iK '+y)~^ 


Moreover, since 

Z{y)=E{y) — y ^ 


there results 

1 1 {K+ iK',y) - II (K, y) = iK'E(y) 

and 

h-=—K + jJ I (98) 

The expression for equation (96) can now be written 



1 2[k- ['i'W -r )[] I -»WA” 

(99) 

where the moduli are k for the non primed functions and 
k' = ^|\ —k~ for the primed functions. By defijiition, 7 = arc 

672, ^ = 7 '’ 7^ where F is the incomplete elliptic integral of 


the first kind with argument 7- and modulus k. 


In the same notation, the equation for Ho is as follows: 


II 0-- 




b^k 

, EK . 
i^{y)—y + 


Ha->)] 


(100) 


Formulas (99) and (100) can now be combined with ecpia- 
tions (79) to give 




and 


„ s ! 2G 


( 101 ) 


whei’e 





2G 

^0+^1 


p 1 “h ^0^1 V ( 1 ^0^) (1 ^1^) 

00+^1 


(102) 


(103) 


and K' is the complete elliptic integral of the second kind 
with modulus — 

From equations (87), (101), and (102) 


Pi-Pu ^Ap^2C{d)^ 2 z\d + B 

f T/ 2 ^ (^0 ^) 

P 0 ^ 0 


2a / 2G “ r 
pE'^jeo+^,l 


(do — ) 9-]- 2^o ^i~l 


(104) 


It should be remaiked that the slope of the loading curve 
is zero at ^=0. Figures 19 (a), 19 (b), and 19 (c) show the 


variation of - with j3 tan 6 for values of /3 tan 61 equal to 
a ([ 

0.3, 0.6, and 0.9, respectively, and for 13 tan 5 q equal to 0, 
0.3, 0.6, and 0.9. 
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From equations (49) and (104) the lift eoefficient for a 
type 3 plan form is obtainable. Two cases will be developed 
here: first, when the trailing edge of the wing is perpendicular 
to the stream direction; second, when the trailing edge of 
the wing is perpendicular to the line of symmetry. The 
first configuration may be referred to as a skewed wing while 
the second configuration may be referred to as a symmetrical 
delta wing at an angle of sideslip. Thus for a skewed wing 

Cr,=^-^^{e, + e,)2G (105) 

where G is given by equation (103) and /i*' has the modulus 
— This result agrees with that given by K. C. 
Koberts in an abstract in reference 13. 

For the more practical case of the delta wing at an angle 
of sideslip, figure 12 (c), the lift coefficient can be expressed as 


where A is the angle of sideslip and 2A the angle between the 
leading edges, and G is expressed in terms of aiul di which 
are, in turn, expressed in terms of A and A by the following 
equations 


6y=^ tan (A -A) 1 


(107) 


vSince the pressure distribution has been computed only for 
wings with leading edges beliind the Mach cone springing 
from the apex and with a trailing edge ahead of the Mach 
cones from the wing tips, formula (106) is valid only for 
cases where 

/. + A<90^\ 

A + A<^/i > (108) 

A-A>0 ) 


Cl =*^7 COS A 




tan A 


(106) 


These restrictions are practically always met, however, for 
angles of sideslip likely to be encountered in flight. 
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Figure 20. — Variation of rediiml lift-curve slope with angle of sideslip A for some 


l)lan forms of type .‘1. 


dC^ . 


Equation (106) is plott(‘d in figure 20 where /3 is 


shown as a function of sideslip and A. The figure shows 

dC 

that up to 15° of sideslip ^ remains practically constant. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., April 191^.7. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

TiOngitndiriAl 

X 

X 

Rolling 

L 

Y >Z 

Roll 

d> 

u 

p 

T iSteral 

Y 

Y 

Pitching 

M 

Z >X 

Pitch. 

Vr 

Q 

u 

q 

Normal— 

Z 

Z 

Yawing 

N 

X >Y 

Yaw 


w 

r 


Absolute coefficients of moment 

p ^ ft ^ 

(rolling) (pitching) 


(yawing) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

P 

Diameter 
Geometric pitch 

P 

P 

Power, absolute coefficient Cp= — orTn 
' pn^D^ 

PID 

V' 

Pitch ratio 
Inflow velocity 

a 

Speed-power coefficient 

V, 

Slipstream velocity 

V 

Efficiency 

T 

Q 

T 

Thrust, absolute coefficient 
Torque, absolute coefficient ^5 

n 

Revolutions per second, rps 
Effective helix angle =tan“^^ 2 ~^^ 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower =0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


5. NUMERICAL RELATIONS 

1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 


